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We extend the worldline instanton technique to compute the vacuum pair pro- 
duction rate for spatially inhomogeneous electric background fields, with the spatial 
inhomogeneity being genuinely two or three dimensional, both for the magnitude 
and direction of the electric field. Other techniques, such as WKB, have not been 
applied to such higher dimensional problems. Our method exploits the instanton 
dominance of the worldline path integral expression for the effective action. 
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I. INTRODUCTION 

A remarkable prediction of quantum electrodynamics (QED) is that in the presence of an 
electric field the polarization of the vacuum can lead to the production of electron-positron 
pairs from vacuum. This was initially predicted and estimated for approximately uniform 
as not yet been observed experimentally due to the extreme scales required 
In achieving such extreme field strengths, the inhomogeneity of the field becomes 
important. Much work has been done to estimate the pair production rate for electric fields 
whose direction is fixed and whose magnitude varies in one dimension, either spatial o, fif or 

nnn 

temporal [3, la I2( • These approaches are semiclassical, essentially WKB or its variants. More 
recently, a Monte Carlo worldline loop method has been developed and applied to the vacuum 
pair production problem . In principle, this Monte Carlo method can be applied to very 
general electric fields, but it has so far only been applied to the case of one-dimensional 
spatial inhomogeneities. The worldline instanton method, in which the Monte Carlo sum is 
effectively dominated by a single instanton loo p, w as introduced for constant fields in 
and extended to inhomogeneous fields in [12I Il3^. For one-dimensional inhomogeneities, 
the agreement between WKB methods, worldline instantons and the Monte Carlo results is 
excellent Q. In this paper, we app.y the worldline instanton ,„ethod to electric he.ds with 
multidimensional inhomogeneities. Specifically, we compute the pair production rate for a 
class of spatially inhomogeneous electric fields, in which both the magnitude and direction 
of the electric field vary in two or three dimensional space. It is not clear how to compute 
the pair production rate for such fields using WKB. 

The vacuum pair production rate can be deduced from the imaginary part of the effective 
action T Q]: 

P pr oduction = l-e- 2Imr ^2Imr . (1.1) 

The technical problem is thus to compute Im V for a given background electric field. This is 
nontrivial as the answer is non-perturbative in the field. For example, for a constant electric 
field the leading weak field result (for scalar QED) is 

Iml e 2 E 2 , s 

1/4 loir 6 

For inhomogeneous fields, the simplest approximation is the "locally constant field" (LCF) 
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approximation, in which one replaces the constant E in (| 1.2(1 by E{x), and integrates: 

9 r 2 

Imr LCF ~ — - / d A xE 2 {x)e *m . (1.3) 

167T d J 

In this paper we compute improved approximations to Im T using the worldline instanton 
method. In SecJnJ we review the worldline instanton method for computing the prefactor; 
in Sec. Illl| we present explicit computations of Im T for two-dimensional spatially inhomo- 
geneous electric field backgrounds, and in Sec. IIV| we present explicit computations of Im T 
for three-dimensional spatially inhomogeneous electric field backgrounds. We conclude with 
some brief comments on possible extensions. 



II. WORLDINE INSTANTON METHOD 



The worldline formalism for QED provides a nonperturbative approach to computing T, and 
in particular its imaginary part. The worldline formalism expresses the one- loop effective 



a quantum mechanical path integral 
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(2.1) 



dV 0) J Vxexp 

x(T)=x(0)=x(°) 

Here the functional integral J T>x is over all closed Euclidean spacetime paths x^t) that 
are periodic (with period T) in the proper-time parameter r. These closed paths are based 



at the marked point x^\ whose location is integrated over 20, 21 1. We use the path integral 
normalization conventions of Q|. The effective action T[A] is a functional of the classical 
background field A^x), which is a given function of the space-time coordinates. 
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that for certain 



The worldline instanton method is based on the observation 
classical background fields the quantum mechanical path integral in (j2.1|) may be dominated 
by an instanton configuration, which is a closed path x M (r) satisfying the classical Euclidean 
equations of motion 



x 2^6 F^yXX) x v , 



(/1, !/ = !.. .4) 



(2.2) 



where = d fl A u — d u A fl is the background field strengt 
solutions to (J2.2)) . For uniform jll| and inhomogeneous 



1. Worldline instantons are periodic 
"1 

12j background electric fields, these 
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classical worldline instantons straightforwardly determine the nonperturbative exponential 
factor in ImT. 

However, this exponential factor is only part of the story. For inhomogeneous fields, 
the prefactor is also important. For example, with a spatially inhomogeneous electric field, 
there is a cutoff in the scale of the inhomogeneity, beyond which the pair production rate 
vanishes. Interestingly, this cutoff is reflected in the prefactor and not in the exponential 
factor [^, . Physically, this cutoff arises when the range of the electric field is so small 

that a virtual pair accelerated by the field does not acquire sufficient energy to become real 
particles. 

To compute the prefactor one must compute the effect of the fluctuations about the 
classical worldline instanton path. The strategy for computing the prefactor was explained 
in The basic idea is that the quantum mechanical path integral in (|2.1|) may be 



approximated as 
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Vx exp 

x{T)=x(0)=x(°) 



x 2 M I e" 5 ^']^) 



dr [ — + ieA ■ x 



(4vr) 2 v^tA 



(2.3) 



where A is the operator describing quadratic fluctuations about the worldline instanton path: 

Id 2 d 
A^ = --5^ — + ie F^ u (x) — + ie d^F^x) x p . (2.4) 
2 ar z ar 



An important technical observation |22j, |23j is that Det A can be expressed as the determi- 
nant of a finite dimensional matrix, whose entries consist of certain solutions to the Jacobi 
equations, 

A^Vu = , (2.5) 

evaluated at r = T. Specifically, we define the four independent solutions rf v > to (|2.5jh with 
initial value boundary conditions 

< } (0) = ; r)H(0)=V . (2.6) 

Then Q Q 

Det(A) = detfoW(D] • (2-7) 

This provides a simple and general numerical method for computing the semiclassical ap- 
proximation in ()2.3|) . for a given proper-time T. The next step [l3| is to evaluate the 
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remaining T integral in (|2.1|) using steepest descents at a particular critical value T c . This 
procedure was implemented explicitly in for inhomogeneous fields varying in one di- 
mension, and the agreement with other computational methods was shown to be excellent. 
Here, in this paper, we extend this computational procedure explicitly to higher dimensional 
inhomogeneities. 

A. Numerical method for finding worldline instanton loops 

To be specific, we consider spatially inhomogeneous electric fields, E(x), for which we write 
the Euclidean gauge field as 

A 4 (x) = -ijf(kx) . (2.8) 
Here k is a parameter related to the scale of the variation of the function /. By analogy with 



Keldysh's adiabaticity parameter 



24j for time dependent ionization problems, we define the 



dimensionless inhomogeneity parameter 

mk 

= < 2 ' 9) 

More general fields are characterized by many scale parameters, but the form in (|2.8j) is 
sufficient to illustrate our results. 

The first numerical step is to find the closed worldline instanton path, as a periodic 
solution to the classical Euclidean equations of motions 

Xj = 2eEdjfx 4 , (j = 1,2,3) , (2.10) 



x 4 = -2eEdjf±j . (2.11) 

To solve numerically these equations, we need to start from a point on the solution curve. 
In principle, this step can be done by a tedious search. However, for the cases considered in 
this paper, all the solutions pass through the local maximum of the electric field. Without 
loss of generality, we choose this point to be the origin: x = (0, 0, 0). We exploit the gauge 
freedom to choose /(0) = 0. 

Note that the i± equation (|2.11|) can be integrated immediately to yield 

2eE 

x 4 = ^f(kx) , (2.12) 
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where we have chosen the integration constant to vanish in order to have a periodic solution. 
Now, recall that (|2.2|) implies that x 2 is constant on a classical solution: 

xl = a 2 . (2.13) 

Thus, the magnitude, |x(0)|, of the "3- velocity" is equal to this constant a. But we do not 
know the initial direction of the 3-velocity. We use as shooting parameters the direction of 
the initial 3-velocity. For a one-dimensional problem where / is a function of just one spatial 
coordinate, no shooting is required. For a two-dimensional problem where / is a function 
of two spatial coordinates, one shooting parameter is required. For a three-dimensional 
problem where / is a function of all three spatial coordinates, two shooting parameters are 
required. Given a value of the parameter a in ()2.13|) . we integrate the equations (J2.1()j) and 
(I2.11j) . adjusting the initial direction of the 3-velocity until a closed loop is found. The total 
proper-time T required to close the loop is clearly a function of the parameter a, so 

T — T(a) or a = a(T) . (2.14) 

For an electric field depending on two spatial coordinates, say X\ and x 2 , the loop has 
X3 = 0, and so depends in general on x±, X2 and £4 only. An example of such a numerical 
worldline loop, for f(kx 1 ,kx 2 ) = i^kxi)^+io^x 2 ) 2 ' * s shown i n Fig- d An example of such 
a numerical worldline loop for an electric field depending on all three spatial coordinates, 
with f(kx 1 ,kx 2 ,kx 3 ) = [1+( ^ i) 2^ 1 2 |^ 2 2 ) t+ 1 3 0(fca; 3 ) 2 ] ^ is shown in Fig. H This loop depends on 
all four space-time coordinates, so we plot appropriate cross-sections. Given these worldline 
loops, the pair production rate can be computed numerically as outlined above and in kJ. 
In the following sections we illustrate this in some specific examples. 



III. TWO-DIMENSIONAL ELECTRIC FIELDS 



In this Section we illustrate our worldline instanton procedure in a class of models where 
the electric field is static and only depends on two spatial coordinates. For the 2D space- 
dependent fields of the form (|2.8j) . the fluctuation operator (|2.4j) can be restricted to its 
components in the (xi,x 2 ,x 4 ) plane : 
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(3.1) 
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FIG. 1: Numerical worldline instanton loop for the case of f{kx\ : kx2) = i + (k x ^z+iQ(kx 2 )Z ■ The 
first plot is the 3D parametrized plot of {xi(r), X2(t), Xi(r)}; the second plot is xi(t) vs. X2(t); the 
third plot is x\{t) vs. x^(t); and the last plot is X2(t) vs. x^{t). The parameters used to generate 
these plots are E = 1, k = \/2, a = 0.593, and xi(0) = 0.562408 091043 522 23. 
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FIG. 2: Numerical worldline instanton loop for the case of f(kxi,kx2,kxs) = 

[1+(to0 H 2 |5)t+To( fc , 3 )^ - The first P lot is the 3D Parametrized plot of {*i(t), x 2 (t), x 3 (t)}; the 
second plot is {x\(t), X2(t), £4(7")}; the third plot is {x2(t), x%{t), X4(r)}; and the last plot is 
{^(t), xi(t), X4(t)}. The parameters used to generate these plots are E = 1, k = 1, £i(0) = 0.6, 
x 2 (0) = 0.453 000 684, and x 3 (0) = 0.20173312728. 
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where 

am an 

To compute the fluctuation determinant we use the semiclassical quantum mechanical path 
integral result (|2.7|) . Thus, we need to find solutions to the Jacobi equations A rj = in (j2.5|) . 
satisfying the initial value conditions (j2.fij) . In general, this needs to be done numerically. 
However, there exists a class of fields for which the entire computation can be done semi- 
analytically, providing more insight. Consider fields (|2.8|) with the symmetry 

f(kx 1 ,kx 2 ) = f(kx 2 ,kx 1 ) . (3.3) 

For example, the two-dimensional symmetric electric field example with f(kxi, kx 2 ) = 
(kxi + kx 2 )e~^ kxi ^ 2 ~^ kX2 ^ 2 , is shown in Fig. El As can be seen from the Figure, such an 
electric field is genuinely two dimensional in its inhomogeneity, both for its direction and 
its magnitude. Physically, we choose the form of / so that the electric field is spatially 
localized. Interestingly, for such fields the worldline instanton loop, which is a closed loop 
in three di m en 81 on 8 fc, „, and x 4 ], is plana, This implies that much of the analysis can 
be reduced to the one dimensional case studied in [13J. Note, however, that even though the 
classical loop is essentially 2D, the fluctuations extend non-trivially into the third direction. 

A. Determinant of the Fluctuation Operator 

We write the solution of the Jacobi equation A</> = as 



(3.4) 



Numerically, one simply integrates the Jacobi equations with the initial value conditions 
(12. 6|) . However, we can partially construct these solutions analytically, as follows. There are 
six linearly independent solutions, from which we need to construct the three independent 
solutions satisfying the initial conditions (j2.fi)) . To find all six linearly independent solutions, 
we use two different Ansatze: 

1. 0i = 02- In this case, the Jacobi equation reduces to 

-\h + eEk [f^ + /(Li)] ^ + eEf d,o)^ = , 

<j> 4 = -4eEfW<j>i + v 4 , (3.5) 
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FIG. 3: The electric field E(xi,x 2 ) for the case of f(k Xl ,kx 2 ) = (feci + kx 2 )e-( kxi ^- ( - kx2 ^ . The 
first plot is the vector plot of E as a function of (kx±, kx 2 ); the second plot is the 3D plot of the 
magnitude of E as a function of (kx\, kx 2 ). To generate these plots, we use the parameter E = 1. 
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where v± is an integration constant. There are four independent solutions of this type. 
As in the ID case considered in Q], we can write these four solutions as : 
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xf(t) 
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X?(t) 



J [xf(ttf J [±f( 



x f(t)Y 
2 



(3.6) 



2. 0i = —02, and 04 = 0. In this case, the Jacobi equation reduces to a single second 
order differential equation: 



(3.7) 



There are two independent solutions of this type. These are new to the two-dimensional 
case. Denote as ip^ and ifj^ 6 ' the two linearly independent solutions of ()3.7|) with the 
initial conditions: 



^ 5) (0) = 1 , 

^ (6) (o) = o , 



^ (5) (0) = 
^ (6) (0) = 1 



(3.8) 



Then, we can write the last two solutions of the Jacobi solutions A0 = as 

/ ^(5) \ 

, (6) (r) 



(5) (r) 



_^(5) 

V / 



/ ^,(6) \ 



_^,(6) 

V / 



(3.9) 
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Finally, given these six linearly independent solutions, cffi^ , . . . , <jfi\ we construct the 
linear combinations satisfying the initial conditions (|2.fij) as 



ixf(O)0( 3 )(r) + i^(r) 



(3.10) 



A lengthy but straightforward computation shows that the fluctuation determinant (|2.7J1 
takes the following simple form: 



Det(A) = Tdet [r/ (1) (T), r/ 2) (T), r, {i) (T)] 



T[x?(0)} : 



■^\T) [a 2 I 2 x {T) - li(T)] 



(3.11) 



2 xf(T) 

Here, I\(T) and Iz(T) are integrals defined in terms of the closed worldline instanton loop: 
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(3.12) 



221. 



The overall factor of T is from the free X3 direction 

The determinant (|3.11|) can be simplified further using periodicity, which implies that 
xf(0) = xf(T), and the vanishing of hiT). Therefore, we find the simple expression for the 
fluctuation determinant: 

2 



Det(A) = -^\T) 



2eE 



x?(0)j(T)h(T) 



Here we have defined 



(3.13) 



(3.14) 



Now recall from (|2.1|) and ()2.3|) that we still need to evaluate the 4-dimensional space-time 
integral over the fixed point on the closed loops: 



dxi(0) dx2(0) dxs(0) g?£ 4 (0) 
(LT) I dx 1 {0)dx 2 {0) 



(LT)j I drxf(0) 



(3.15) 
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where L is the 1-space volume and T is the total time. N is a normalization factor that 
depends on the orientation of the loop. The simplest way to evaluate N is to compare the 
k — ► limit with the corresponding locally constant field approximation answer (jl.3)l . as 
explained below. Observe that the factor of xf(0) appearing in (J3.15)) cancels against the 
same factor in ^/Det(A) in (j3.13j) . so that the spacetime integration effectively contributes 
a volume factor (LT), a factor of N/k, and a factor of -| from the r integral. This last 
factor is just the collective coordinate contribution arising from invariance under shifts of 
the starting point on the loop, which gives rise to the second of the zero modes in (|3.6j) . as 
discussed in the ID case in 13j. Thus, collecting all the pieces, we see that the worldline 
instanton approximation (J2.3|) to the quantum mechanical path integral leads to : 

Y^^LtJL VjL r dT exp [~ {S[±*](T) + m?T)] 
(4vr) 2 AeE J y/r j(T)h(T)^M(fj 

The main difference from the ID case in [13j is the appearance of the factor of l/y/^{T) 
and the power of T. The former is a measurement of the fluctuation in the second spatial 
direction; and the latter is because the remaining free dimension is now one instead of two. 

B. The T integral 

As in [3], for spatially inhomogeneous electric fields we use a rotated steepest descents 
method to evaluate the T integral in the vicinity of a critical point T c . The critical point is 
a stationary point of the exponent in ()3.16|) : 

A(T) = S[x cl ](T) + m 2 T . (3.17) 

This notation emphasizes the fact that the action 5"[x cl ], evaluated on the worldline instanton 
path x cl (r), is a function of T. Near T c , we have 

A(T) « A(T C ) + ±A"(T C )(T - T c ) 2 + ■ ■ ■ (3.18) 

As in the ID cases, the critical point T c occurs when 

7(T C )=7 • (3.19) 

The rotation of the T contour produces a factor of i, so our final expression is 

N 1 P -A(T C ) 
.wli . . f L7 -\ iV _± e C3 20) 

32tt 3 / 2 eE ^(TJV^WIA"^)! 
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C. Locally Constant Field Approximation 

As in jfj we compare the worldline instanton result with the locally constant field 

approximation (jl.HJ) . For example, consider the explicit example 

kxi + kx 2 



f(kx 1 ,kx 2 ) 



(3-21) 



1 + (k Xl ) 2 + (kx 2 ) 2 ' 

Expanding \E(xi, x 2 ) | for small k, we find it is a quadratic function of a new pair of variables: 



where x + = X\ + x 2 , and x 

imr LCF 



/ 3 

E{xi,x 2 )\ ^V2E[l--k 2 x 2 + --k 2 xl + -- 

\ x\ — x 2 . Thus, the leading LCF approximation is 



(3.22) 



(LT) ^ o / dx + dx_ exp 



32tt 3 



m 2 7r 



V^eE V 2 



-k 2 x 2 



(LT)- 



[eEf 



-g V5e£ 



(3.23) 



4V37r 3 (mk) 2 

Comparing Imr LCF and the limiting value of Imr WLI as k — > 0, we find empirically that the 
normalization constant in ()3.2()j) for the case of ()3.21|) is 

r eE 



N = 2 5 / 4 



in 



(3.24) 



The same normalization factor arises for the exponentially localized field with f(kxi, kx 2 ) = 
(kx\ + kx 2 )e~( kxi ^ 2 ~( kX2 ^ 2 . The ratio of the worldline instanton answer (J3.20j) to the LCF 
answer ()3.23|) is plotted in Fig. 0] for two different forms of the 2D electric field. In each case, 
there is a smooth dependence of the ratio from 1 to as the inhomogeneity scale 7 [i.e. k] 
increases, which corresponds to the range of the field decreasing. Physically, the vanishing of 
the imaginary part at a critical inhomogeneity scale is because for very short range electric 
fields virtual dipole pairs do not acquire sufficient energy when accelerated in the field in 



order to become rea 
inhomogeneities 



particles. This has been observed previously for one dimensional spatial 
I, and here we see it also in the higher dimensional case. 



IV. THREE-DIMENSIONAL ELECTRIC FIELDS 

For three dimensions the computation is very similar, with just a few changes. First, finding 
the worldline instanton loop requires a two parameter shooting procedure. Second, the 
fluctuation operator A in ()2.4|) is now a 4 x 4 matrix differential operator. The computation 
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Q- -El f(kx ,kx )=(kx +kx )e 



G-0 f(ix 1 ,kx 2 )=(kx 1 +kx 2 )/[l+(kx 1 ) +(kx 2 )")] 




FIG. 4: Imr WLI /Imr LCF for the cases of f(kx u kx 2 ) = l+[ ^$ k + x { kx2) ^ and f(kx u kx 2 ) = (kx\ + 
kx2)e~( kxi } 2 ~( kx2 ^ 2 . The parameters used to generate this plot are E = 1, e = 1, and m = 1. 

can be done numerically, but as in the 2D case we can gain more insight by specializing to a 
class of fields for which Im T can be computed semi-analytically. Consider totally symmetric 
fields (£21) with 



f(kxi, kx 2 , kxs) = f(kx2, kxi, kxs) = f(kxi, kx%, kx 2 ) = f(kx-s, kx 2 , kx±] 



(4.1) 



v = is shown in 



The electric field for an example with f(kx±, kx 2 , kx 
Fig. El Note that it is localized, and that both the magnitude and direction of the field vary 
in three dimensions. 

The calculation is similar to the two-dimensional symmetric case, except we now need 
eight [rather than six] independent solutions to the Jacobi equation, A0 = 0, in order 
to construct the four [rather than three] independent solutions 77 ^ satisfying the initial 
value conditions (|2.6[) . In order to find the extra two solutions we extend the Ansdtze for 

= (01,02,03,04) to: 



1. Suppose 0i 



'3- 



Then A0 = reduces to: 



-|0! + eEk [f^ + 2/( 1 > 1 >°)] xffa + eEkf 1 



mi 



U = -6eEf 



(1.0,0), 



(4.2) 
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FIG. 5: Vector plot of the electric field E(x±, X2, X3) for the case of f(kxi,kx2,kxs) = 

kxi+kx2+kx3 

[l + (fcXl)H(fc2 2 ) 2 + (fcx 3 ) 2 ] 3/2 ' 

There are four solutions of this type. These solutions of the Jacobi equation are similar 
to those in the one- and two-dimensional cases. 

2. Suppose 4>\ = —</>2 and 03 = <p4 = 0. In this case, the Jacobi equation reduces to a 
single second order differential equation, which is very similar to (J3.7)) : 

-& + eEk[f<WM-f<- 1 > 1 M]£&fa = o . (4.3) 

There are two solutions of this type. Define ip^ and to be the solutions to (|4.3|) 
with the initial conditions in (|3.8j) . Then two further solutions can be written as 



( ^ (5) \ 






(6) (r) 



I ^ (6) \ 

_^(6) 







(4.4) 
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3. Suppose 0i = 02 = — |03j and 4 = 0. In this case, the Jacobi equation reduces to 
exactly the same equation as in ()4.3|) . There are two further solutions of this type. 
They can also be written in terms of ip^ and if)^: 

( ^ \ ( ^ (6) \ 



0( y )(i 



V 



^(5) 





r) 



^(6) 





(4.5) 



Given these eight linearly independent solutions to A0 - 
satisfying the initial conditions ()2.6jl are 

|xf(O)0( 3 )(r) + |0( 6 )(r) + | 



0, the linear combinations 



V {1 \r) 



V i2) (r) 



ixf(O)0( 3 )(r)-i0( 6 )(r) + |0( 8 ) 
Ixf(O)0( 3 )(r)-i0( 8 )(r) , 



|^(O)0( 3 )(i 



3 r ■ (4-6) 

A lengthy but straightforward computation shows that the fluctuation determinant ()2.7j) 
takes the simple form: 



Det(A) = det [r/ (1) (T), r/ 2) (T), r/ 3) (T), r/ 4) (T)] 



(4.7) 



where I\ and J2 are defined in ()3.12|) . Notice the remarkable similarity to the ID and 2D 
results for Det (A), in [ljj and ()3.13j) . respectively. 

Proceeding as in the 2D case, we find the following worldline instanton expression for the 
imaginary part of the effective action: 

N 1 e~ A ^ 



ImT 



WLI 



rJl 



2 32^k eEjl^TWHTc) V\&"<M ' ^ 
Here A^ is a normalization constant coming from the spatial integrals, which depends 
on the form of the planar worldline instanton loop. As before, the simplest way to 
evaluate A" is by comparison of the k — > limit of Im T WLI with the locally constant 
field approximation expression in (jl.3j) . For example, for the case of f(kxi, kx 2 , kx%) = 
(kxi + kx 2 + /c^e - ^ 1 ) 2- ^ 2 ) 2- ^ 3 ) 2 we find 



Im r 



LCF 



^ 3 9 / 4 {eEfl 2 

T - — - e v^ eE 

16\/7 7r 3 (m/c) 3 



(4.9) 
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B -H f(kx r kx 2 ,kx 3 )=(kx 1 +x 2 +x 3 )e" (kx ' ) " (kX2) " <k! ' 3r 
G-O f(kx , ,kx 2 ,kx 3 )=(kx ] +kx 2 +kx 3 )/[ 1 +(kx 1 ) 2 +(kx 2 ) 2 +(kx 3 ) 2 f P ~ 




0.8 



kxi+kx2+kxs 



and 



FIG. 6: Imr WL1 /Imr L " for the cases of f(k Xl ,kx 2 ,kx 3 ) = 3/2 : 
f(kx\,kx2,kx%) = (kx\ + A;x2 + A;x3)e _ ^ fca:i ) 2_ ^ c:r2 ) 2_ ^ :r3 ^ 2 . The parameters used to generate this 
plot are E = 1, e = 1, and m = 1. 



kx\+kx2-\-kx3 



which leads to N = -4= ^4. And for the case of f(kxi, kx2, kx 3 ) = - — - 
we find 



3 3 / 4 (e£) 



7/ 2 



Imr LCF ~ T— — v ;7 y ..„ e~^ 



4\/l4 Ti 3 (mky 



(4.10) 



and iV = -| 4. 

In Fig. IH1 we plot as a function of the inhomogeneity scale 7, the ratio of the worldline 
instanton answer (|4.8jl to the LCF approximation answer (|4.9jl or (|4.1()jl . for these fields. 
Note once again that the pair production rate vanishes when the field becomes too closely 
localized. 



V. CONCLUSIONS 

In this paper we have extended the worldline instanton method of jl^l . to the multi- 
dimensional case. Specifically, we have computed the imaginary part of the effective action 
in scalar QED for background electric fields that are inhomogeneous, both in their direction 
and magnitude, in two and three spatial dimensions. As in the one dimensional cases, the 
results exhibit the vanishing of the imaginary part when the inhomogeneity scale becomes too 
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small, consistent with the physical picture of pair production as the acceleration of virtual 
dipole pairs in the electric field. In these multi-dimensional cases, there are currently no other 
computations with which to compare, except the crude locally constant field approximation, 
which does not capture the physics of the vanishing of Im T at short inhomogeneity scales. 
It would, however, be very interesting to compare with the numerical Monte Carlo worldline 
loop approach of Gies et al [ll]] , to see if the instanton dominance of the worldline instanton 
method can be combined with the versatility of the Monte Carlo approach. Other interesting 



open questions include: (i) more general Abelian fields representing realistic laser fields 

hn 

with both spatial and temporal inhomogeneities; (ii) incorporating kinetic effects [25, 26] and 



271; 



backreaction effects 
for heavy ion physics |29 
equations 



in) the extension to non- Abelian theories [2j|, such as are relevant 
|, and where the worldline instanton equations are the Wong 
In these more general cases, finding the worldline instanton closed loops will 
most likely require going beyond simple shooting techniques, to use action minimization 
algorithms. 
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